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Abstract 

A dynamical model of the dielectric screening of conjugated polymers is introduced and solved 
using the density matrix renormalization group method. The model consists of a line of quantized 
dipoles interacting with a polymer chain. The polymer is modelled by the Pariser-Parr-Pople (P-P- 
P) model. It is found that: (1) Compared to isolated, unscreened single chains, the screened l l B~ 
exciton binding energy is typically reduced by ca. 1 eV to just over 1 eV. (2) Covalent (magnon 
and bi-magnon) states are very weakly screened compared to ionic (exciton) states. (3) Screening 
of the exciton is closer to the dispersion than solvation limit. 

PACS numbers: 78.67.-n, 77.22.-d 
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I. INTRODUCTION 



The intra-chain excitations of a conjugated polymer in the solid state are often strongly 
screened by the environment of the other polymers. This screening causes significant re- 
ductions of the ionic single-chain excitation energies and charge gaps, and thus single-chain 
calculations are expected to considerably over-estimate the actual excitation energies of 
polymers in the solid state. In particular, the energy of a point charge is typically reduced 
by approximately 1 eV, while the excitation energy of the lowest dipole allowed exciton is 
reduced by approximately 0.5 eV, implying a deviation of approximately 1.5 eV in the single- 
chain calculated exciton binding energy in comparison to the solid state measurements.—^ 

It is useful to consider two limiting cases for the interaction between the electron-hole pair 
and the polarization of the dielectric medium, namely solvation-like versus dispersion- like 
interactions. By solvation, we mean a polarization that develops around the electron and 
hole and follows the motion of these charged species (that is, the dielectric is fast compared 
to the electron-hole motion). In this case, the effects of the polarization can be absorbed into 
a static screened electron-hole Coulomb interaction. By dispersion, we mean the interaction 
between fluctuating dipoles, such as may arise between the fluctuating dipole of the bound 
electron-hole pair in an exciton and the dipoles of the dielectric medium. The interactions 
become dispersion-like when the time scale of the electron-hole motion is comparable to or 
less than the time scale of the dielectric response (that is, the electron-hole motion is fast 
compared to the dielectric). In this case, the dielectric polarization does not identically 
follow the electron-hole motion, but rather fluctuates in a manner that is correlated with 
the motion of the electron-hole pair. Now, the effects of the polarization may be modelled 
by a dynamically screened electron-hole Coulomb interaction. 

The nature of the screening thus depends on the relative time scale of the electron-hole 
motion, which is approximately inversely proportional to the exciton binding energy 3 , versus 
the time scale of the dielectric polarization, which is inversely proportional to the optical 
gap of the dielectric medium. In inorganic semiconductors, the exciton binding energy is 
one to two orders of magnitude smaller than the optical gap. The polarization is then much 
faster than the electron-hole motion, such that it can follow the electron-hole motion and 
effectively screen the Coulomb interaction.— In organic semiconductors, the exciton binding 
energy is comparable to the optical gap and we expect the interaction between the exciton 
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and the dielectric polarization to have a dispersion- like character. In this regime, there is no 
a priori reason to suppose that a static dielectric response function will correctly describe 
the relevant physics 2 . 

In this paper we present a dynamic model of the dielectric - making no assumptions about 
time scales - with which we attempt to understand the screening of the key excited states 
of conjugated polymers. We consider the role of lengths scales (namely, the electron-hole 
separation versus the length scales in the monopole-dipole interaction), time scales (namely, 
exciton excitation energy versus binding energy), and ionicity versus covalency in the excited 
states. The correlated electron dynamics within a single chain are modelled by the Pariser- 
Parr-Pople (P-P-P) model, while the dielectric is treated as a linear chain of quantized 
dipoles. This model is solved essentially exactly via the density matrix renormalization 
group (DMRG) method.- 

The plan of this paper is as follows: in the next section the model is introduced and 
parameterized. In the following section we discuss our results, and finally conclude. 

II. DYNAMIC DIELECTRIC MODEL 

The model is designed to represent the dielectric response of an array of conjugated 
polymers parallel to the solute (or test) chain. The dielectric function, e(a>), of an assembly 
of polymers is£, 



where u n and f n are the frequency and oscillator strength for the nth transition, q and m 
are the electronic charge and mass, eo is the dielectric constant of free-space, and N is the 
number of polymers per unit volume. In this paper we consider a 'two-state' model, namely 
we only consider excitations of the dielectric from the ground state to the first dipole-allowed 
state (e.g. the l 1 !^ state in centro-symmetric polymers). Then, using the oscillator sum 
rule, 
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FIG. 1: The polymer chain (dashed line), with a monopole at Xi (circle) and the line of longitudinal 
dipoles (dotted line), with a dipole at Xj (arrow). 



we have 

e((jj) = 1 + — —— 2 rr. (4) 

e h{u 2 - w 2 ) 

In this expression, (/}) and co>o are the transition dipole moment and frequency of the excited 
state. Tiujq corresponds to the optical gap of the polymers. We note that for to » lu Eq. 
(HJ) becomes, 

= 1-4' ( 5 ) 
where u 2 = Nq 2 je^vn is the plasma frequency. This frequency is typically 10 — 20 eVs - 
much larger than the optical gap or the energy scale of the electron-hole motion - so this 
limit need not concern us here. 

Since the dipole moment of the chain is parallel to the chain axis, we consider a line of 
longitudinal dipoles, each oscillating with a characteristic frequency, c^o-" The dipole at Xj 
interacts with a point-charge on the solute chain at x« via the monopole-dipole interaction, 

Us{X ' r) = Ane (xZr^ (6) 
where /i is the dipole moment, ) and r is the normal distance of the line of 

dipoles from the chain, as illustrated in Fig. 1. 

A realistic model of the dielectric would consist of a three-dimensional array of lines of 
dipoles surrounding the solute chain. However, such a model would not be easily soluble by 
the DMRG method. Instead, we consider a one-dimensional model and derive an effective 
interaction between the monopoles and the three-dimensional array of dipoles by integrating 
over a cylinder with inner and outer radii of r\ and r2, respectively: 

U(x; r 1; r 2 ) = / p2narU s (x, r)dr (7) 
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FIG. 2: The integrated monopole-dipole interaction, U(x;ri,r2), for n = 4 A, with r<i = 8 A 
(squares) and r 2 = 16 A (circles). 

= (dU (tt) H^) h ( tan_1 G9) - sin ( tan " (£)). • 

Here, p is the number density of dipoles and a is the linear separation (for convenience 
taken to be the separation between the sites on the chain). The purpose of this integration 
procedure is that it performs a type of average over the distances between the line dipoles 
and the solute chain. However, r\ is still the effective length scale in the interaction, as 
shown in Fig. 2, which shows U(x; r 1 ,r 2 ) for different r 2 for fixed r±. For both r 2 = 8 A and 
16 A U{x;r\,r2) is strongly peaked at approximately the value of r\. Thus, using a single 
line dipole a distance r\ away would give very similar results to the cylinder of inner radius 

n. 

The Hamiltonian is thus, 

H = H e + H d + H e _ d , (8) 
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where, 

H e = - £ U(clci+ia + 4 +1(rCia ) + fa - I) (Nn - ±) + \ £ ViA^i - l)(i% - 1)(9) 

is the P-P-P model, describing the interacting electronic degrees of freedom. c\ a creates an 
electron with spin a in the 7r-orbital on site i and iVj = J2a c \a c i<r- Vij is the Ohno-Coulomb 
repulsion, 

v a = u = n , . ( 10 ) 



where e is the dielectric constant (generally taken to be 1), the bond lengths are in A and 
j3 = (?7/14.397) 2 . The bond lengths, a, used in the evaluation of are 1.4 A, and the bond 
angle is 180°. The transfer integral is £j = £(1 + (—l) l 5), where 5 is the bond dimerization 
parameter. 

The Hamiltonian, 

H d = Y,^ b %, ( n ) 

describes the oscillating dipoles. fej is a boson operator creating a quantum of energy Tiujq in a 
linear harmonic oscillator located at Xj. This term describes the non- interacting self-energy 
of the oscillating dipoles. However, in this model we neglect the dipole-dipole interactions, as 
in one dimension this induces an unphysical spontaneous dipole moment. Thus, the classical 
limit (defined as uo — > 0) is not completely described by this model. 
Using Eq. (J7J), the monopole-dipole Hamiltonian is, 



H e _ d = 



q2 ^ ^ r._^._-i^^ .-^- 1 /^m ( ^_ 1)a M±M ; 



sin j tan 1 ( — jj — sin ^tan 1 y — 



(12) 



47re / \ r\ 

where we define the dimensionless monopole-dipole coupling constant, A, as, 

X = prl(d) (13) 

and we have used 

A = 5 <<i>M±M (14) 

q(d) is the transition dipole moment of the solvent chains associated with the exciton^ 
and (bj + bj)/2 is the dimensionless operator corresponding to the displacement of the os- 
cillator. 
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In the absence of electron-electron interactions and for only local monopole-dipole inter- 
actions the Hamiltonian, H, is the Holstein model, which is widely used to study charge- 
density-wave phenomena in molecular solids. The monopole-dipole interactions imply that 
H possesses a special type of particle-hole symmetry. In particular, H is invariant under 
the simultaneous particle-hole transformation o& 



hu is the excitation energy of the 1 X B U exciton and q(d) is its corresponding transition 
dipole moment. These are determined by solving the single chain P-P-P model in the absence 
of monopole-dipole interactions, p is the only adjustable parameter in the model, and is 
adjusted to reproduce the solvation energy of a point charge found by Yaron and Moore^ for 
particular values of the P-P-P parameters, namely, t, U and 5. Having fixed p for one set 
of P-P-P parameters the model can be transferred to another set of P-P-P parameters, by 
recalculating the unscreened values of Huq and (d) , and the solvation of states as a function 
of these parameters can be monitored. 

A. Parametrization of the model 

Moore and Yaron^ studied solvation effects in polyacetylene by surrounding a central 
solute chain by increasing numbers of solvent chains, arranged as in the crystal structure 
of refi^, and extrapolating to an infinite system. The central chain was treated with a P- 
P-P model with parameters, t = 2.4045 eV, 5 = 0.0734, U = 11.13 eV (corresponding to 
td = 2.581 eV and t s = 2.228 eV) and the solvent chains were treated in an independent- 
electron approximation. Coulomb interactions were included between all chains. Replacing 
the central solute chain with a point charge led to an estimate of about 1 eV for the solva- 
tion energy of a point charge. In the assumption of an infinitely fast dielectric, where the 
solvent polarization is equilibrated to the instantaneous position of the electron and hole, 
the solvation energy of a well separated electron and hole is about 1.9 eV. Using a model 
that approximates the dynamic response of the dielectric, the free electron and hole become 




(15) 



and a reversal of parity in the boson operator, 



b] - -b). 



(16) 
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TABLE I: Unscreened energies, the transition moment, (d), and the dipole 'polarizability', a 

Id) 2 

of a 50-site polymer chain modelled by the P-P-P model. The energies are in eVs. 



E{V-Bu) 



t U 5 E(1 1 S-) l x B~ binding energy E(2 1 A+) E(l 3 £+) (d) A a 
2.4045 11.13 0.0734 2.85 2.29 2.31 1.28 8.72 26.7 



dressed by the dielectric response of the solvent chains to form polarons and the solvation 
energy drops to about 1.5 eV. Using the ab initio DFT-GWA-BSE method, Bobbert and 
co-workers calculated a point charge screening in polythiophene of 0.55 eV, i.e. 1.1 for the 
charge gapAi. 

Notice that for the parameter values used by Moore and Yaron the P-P-P model is 
strongly correlated^. The triplet, 1 3 -B„ , state is more properly described as a magnon than 
the n = 1 Mott-Wannier triplet exciton. Similarly, the 2 1 A£ state is not the n = 2 Mott- 
Wannier singlet exciton, but a bi-magnon, whose energy lies below the n = 1 singlet exciton, 
namely the l l B~ state . 13 ' 14 The unscreened energies of a 50 site chain are shown in Table I. 

We solve Eq. (jSJ) by the DMRG method 5 using a program well-tested on similar electron- 
boson models 15 . Since the boson energy scale, hu>o, is comparatively large compared to 
the electronic excitation energies, we found that only one boson per site is required for 
the excitation energies to converge to within 0.01 eV. The number of superblock states 
was typically ~ 100, 000 and one finite-lattice sweep at the target chain size was always 
performed to ensure convergence. 

Taking r x = 4 A, r 2 = 16 A and A = 0.14, we obtain a screening energy of 0.81 eV for 
the point charge. This is intermediate between the values of BobbertjAA and Moore- Yaron^ 
and is close to the charged-polaron result of 0.88 for 50 sites.— The distance dependence 
of the charge-charge screening arising from the dielectric model of Eq. ((7j) is similar to 
that obtained from the explicit solvent model of Moore and Yaron 2 , suggesting that the 
longitudinal dipoles provide a reasonable description of the dielectric medium. 

The excitation energy of the l 1 -^ exciton is reduced by 0.69 eV, which is roughly twice 
as large as the Yaron-Moore result^, and much larger than the DFT-GWA-BSE calculation, 
which predicts a negligible screening of the exciton^. The binding energy is therefore reduced 
from 2.3 eV to 1.4 eV. As already noted, for these parameter ranges the l 3 B^ state is a 
magnon and the 2 1 v4+ state is a bi-magnon. Both states are therefore strongly covalent with 
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TABLE II: Screened energies, E, screening energies, S, and the screened 1 B u binding energy, 
BE, of a 50-site polymer chain. The charge gap is denoted by 2A. The energies are in eVs. 



t U 6 A E(l l B-] 


1 S^B-] 




) S{2 1 A+) E(1 3 B+) S(1 3 B+) 2A S(2A) BE 


2.4045 11.13 0.0734 0.14 2.16 


0.69 


2.22 


0.09 1.25 0.03 3.53 1.62 1.37 



very little ionicity and thus there is very little screening. This affect can be traced back to 
the monopole-dipole interaction, Eq. (|12j) . which shows that the electric dipole couples to 
charge density (or ionic) ffuctuations in the wavefunction. These results are summarized in 
Table II. 

III. RESULTS 

Having parameterized the model, we now investigate other parameter ranges. A is found 
using Eq. (JT3J) for fixed p and r\ with the recalculated value of (d). First, we consider 
U = 10 eV, t = 2.5 eV and 5 = 0.1, as this corresponds to the optimal parameter set 
used by Barford and Bursill in their study of trans-polyacetylene^. The single chain P-P-P 
excitation energies are shown in Table III. Notice that, unlike the case for the original P- 
P-P model parametrization, the 2 1 v4+ state lies higher in energy than the l l B~ state. The 
screened energies are shown in Table IV, where we see that the charge gap is solvated by ca. 
2.0 eV, and the 1 X B~ and 2 1 A£ states are screened by 0.80 eV and 0.75 eV, respectively. 
The binding energy of the exciton is reduced from 2.3 to 1.1 eV. 

We next investigate the role played by the strength of the Coulomb interactions, partic- 
ularly on the screening of the 2 1 A+ state, which evolves from a bi-magnon to the n = 2 
Mott-Wannier singlet exciton, and the triplet, 1 3 B^, state, which evolves from a magnon to 
the n = 1 Mott-Wannier triplet exciton, as U is reduce d 13 ' 14 . We take a value for 5 of 0.2 as 
this maps the linear chain band structure onto the band structure of poly(p-phenylene). Ta- 
ble IV confirms the expected trend that as the ionicity of the 2 1 A+ and l 3 B+ states increases 
relative to their covalency their screening, relative to the ionic 1 X -B~ state, increases. 

Length scales also play a role in the screening of the excited states and the charge gap. As 
the separation between a polymer and a line of dipoles, r, is increased the typical distance 
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TABLE III: Unscreened energies, the transition moment, (d), and the dipole 'polarizability' 

(d) 2 

a = )' of a 50-site polymer chain modelled by the P-P-P model. The energies are in eVs. 

E(1 L B U ) 

t U 5 E(1 1 B-) l 1 ^- binding energy E(2 1 A+) E(l 3 5+) (d) A a 
2.5 10 0.1 2.67 2.32 2.95 1.69 9.41 33.2 

2.5 10 0.2 3.69 3.05 4.69 2.76 8.15 18.0 

2.5 8 0.2 3.20 2.76 4.42 2.82 8.81 24.3 

2.5 6 0.2 2.93 2.42 4.10 2.79 9.14 28.5 



TABLE IV: Screened energies, E, screening energies, S, and the screened 1 1 B U binding energy, 
BE, of a 50-site polymer chain. The charge gap is denoted by 2A. The energies are in eVs. 



t U 5 A E(1 1 J B-) S(1 1 J B-) E(2 1 A+) S(2 1 vl+) E(l 3 £+) S(l 3 B+) 2A S(2A) BE 



2.5 10 0.1 0.151 


1.87 


0.80 


2.20 


0.75 


1.59 


0.10 


3.00 


1.99 


1.13 


2.5 10 0.2 0.131 


3.22 


0.47 


4.34 


0.35 


2.73 


0.03 


5.43 


1.31 


2.21 


2.5 8 0.2 0.141 


2.60 


0.60 


3.45 


0.97 


2.71 


0.11 


4.26 


1.70 


1.66 


2.5 6 0.2 0.147 


2.19 


0.74 


2.48 


1.62 


2.40 


0.39 


3.26 


2.09 


1.07 



in the monopole-dipole interaction increases. When this interaction distance is larger than 
the particle-hole separation, the charge character of the exciton approaches that of a point 
dipole. For the line dipole used here, the ratio of the screening energy of the \. l B~ state 
versus the charge gap is 0.6 when r = 4 A and reduces to a roughly constant value of 0.4 
at r > 16 A. This indicates that a substantial fraction of the screening energy of the l 1 -^" 
exciton is from dipole-dipole interactions. 



IV. CONCLUSIONS 



In this paper we introduced a dynamical model of a dielectric to investigate the screening 
of intra-chain excitations by the environment. This model consisted of a linear array of 
quantized longitudinal dipoles interacting with a conjugated polymer, represented by the 
Pariser-Parr-Pople model. The dipole parameters are adjusted to provide a realistic model of 
the dielectric. The model reproduces the expected solvation energy of a point charge in these 
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materials and sets the frequency of the dielectric response to the optical gap of the material. 
The model is one-dimensional and is solved by the density matrix renormalization group 
method. The one-dimensionality of the model does mean that it is impossible to include 
dipole-dipole interactions, and thus the model misses some of the physical characteristics of 
real dielectric materials. 

We make the following conclusions from our results. (1) The 1 X B~ exciton binding 
energy is typically reduced by ca. 1 eV to just over 1 eV. Thus, this relatively simple one- 
dimensional model of the screening of conjugated polymers by the environment predicts 
l l B~ binding energies rather close to the experimentally accepted solid state values of ca. 
1 eV. (2) Covalent (magnon and bi-magnon) states are very weakly screened compared to 
ionic (exciton) states. (3) Screening of the l x B~ exciton is closer to the dispersion than 
solvation limit. 
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